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ABSTRACT 

• i-H ' We take advantage of the exceptional photometric coverage provided by the combi- 

^ \ nation of GALEX data in the UV and the ALHAMBRA survey in the optical and 

near-IR to analyze the physical properties of a sample of 1225 GALEX-selected Ly- 
man break galaxies (LBGs) at 0.8 < z < 1.2 located in the COSMOS field. This is the 
largest sample of LBGs studied at that redshift range so far. According to a spectral 
energy distribution fitting with synthetic stellar population templates, we find that 
LBGs at z ~ 1 are mostly young galaxies with a median value of 341 Myr and have 
intermediate dust attenuation, (E S (B — V)) ~ 0.20. Due to their selection criterion 
they are UV-bright galaxies and have high dust-corrected total SFRs, with a median 
value of 46.4 M©yr _1 . The median value of the stellar mass of the LBGs in the sample 
is log M*/Mq = 9.74. We obtain that the dust-corrected total SFR of LBGs increases 
with stellar mass and the specific SFR is lower for more massive galaxies (downsizing 
scenario). Only 2% of the galaxies selected through the Lyman break criterion have 
an AGN nature. LBGs are mostly located over the blue cloud of the color-magnitude 
diagram of galaxies at their redshift, with only the oldest and/or the dustiest devi- 
ating towards the green valley and red sequence. Morphologically, 69% of LBGs are 
disk-like galaxies, with a fraction of interacting, compact, or irregular systems being 
much lower, below 12%. LBGs have a median effective radius of 2.48 kpc and bigger 
galaxies have higher total SFRs and stellar masses. Comparing to their high-redshift 
analogues, we find that LBGs at lower redshifts are bigger, redder in the UV contin- 
uum, and the presence of older stellar populations in their SEDs is more remarkable, 
although there is no significant difference in the distributions of stellar mass or dust 
attenuation. 
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1 INTRODUCTION 

Much effort has been devoted to searching and analyz- 
ing the physical properties of high-redshift star-forming 
(SF) galaxies. Different selection criteria give rise to the 
discovery of different kinds of galaxies. Among those 
criteria, the most successful and used are the Lyman- 
alpha and the Lyman break techniques, which pick up 
Lyman-alpha emitters (LAEs) and Lyman break galax- 
ies (LBGs), respectively. The Lyman-alpha technique is 
based on looking for a Lyman-alpha emission in the red- 
shifted optical spectrum of galaxies by employing a com- 
bination of narrow and broad-band filters. Specifically, the 
narrow-band filter is used to isolate the Lya line and 
the broad band one(s) to constrain its nearby contin- 



uum dCowie fc HrJl998l;lGronwall et al.l2007l;lGawiser et ail 
20061: lOuchi et al l 120081 . l20ld : iBongiovanni et all l20ld : 
Shiova et all 120091 ). The choice of the central wavelength 



of the narrow-band filter determines the redshift of the 
selected LAEs, which have been se arched / found from 



Nilsson et all 20071. 20091; 


Ouchi et al. 20081; iGawiser et al.l 


20061; Finkelstein et all 


2009b; 


Muravama et al.l 12007; 


Pirzkal et al.l2007; Hibon et al.l201ll. 20121). LBGs are found 



by employing a combination of broad-band filters that sam- 
ple the red-ward and blue-ward zones of the redshifted Ly- 
man break of galax i es, located in 912 A in t he rest-frame 
jMadau et al.l Il996l ; ISteidel et all Il996l . I2003T I. The choice 
of the red-ward and blue-ward filters determines the lo- 
cation in wavelength of the Lyman break and also con- 
strains the redshifts of the selected galaxies. Many samples 
of LBGs have been foun d and examined at different red- 
shifts, mostly at z > 3 dMadau et al.l 19961; Steidel et all 



199d(l999l,[2003l;IStanwav et al.ll2003l;lGiavalisco et alj2004 
Bunker et al.l 12004 IVerma et al.l " l2007l ; llwata et al.l I2007D . 
At lower redshifts (z < 2) the number of LBGs re- 
ported and studied is much lower than that at high e r red- 
shifts llBurgarerla et all 120061. 120071: iLy et all 120091. 1 2011 ; 



Basu-Zvch et al.l201ll;lHathi et al.ll20ld;lNilsson et al.lboil 



Hathi et al]|2012i ; iHaberzettl et al.1 12012 ) despite this red- 
shift range is quite important since it is when it is thought 
that the peak of the cosmic star formation of the universe 
took place. 

Apart from Lya and Lyman break techniques, some 
other methods have been reported to segregate high-redshift 
galaxies in the literature. With the aim of finding galaxies in 
the re dshift gap between z ~ 1 and z ~ 3, lAdelberger et all 
i|2004l ) defined some color selection criteria by employing dif- 
ferent combinations of optical colors to produce samples at 
different redshift ranges: GRi for redshifts 0.85 < z < 1.15, 
GRz for 1.0 < z < 1.5, and U„GR for 1.4 < z < 2.1 and 
1.9 < 2; < 2-7- The galaxies selected in this way have been 
traditionally called BM/BX galaxies. Another ground-based 
optical color selection criterion is the BzK method, which is 
aimed at finding galaxies in the redshift range 1.4 < z < 2.5 
and classifying th em as SF or passively evolving systems 
l|Daddi et a l. 2004) . Both kinds of gala xies are often associ- 
ated to LBGs. IHaberzettl et al.l (|2012T l find that NUV data 
provide greater efficiency for selecting SF galaxies. Further- 
more, they report that, although the BM/MX and BzK 
techniques are very efficient for detecting galaxies within 
1-0 < z < 3.0, the galaxies found are biased against those 



SF galaxies which are more massive and contain a noticeable 
amount of red stellar populations. They argue that, there- 
fore, a NUV-based LBG selection criterion is more adequate 
to compare with populations found at z > 3.0. 

The physical properties of high-redshift SF galaxies 
are usually analyzed by fitting their observed spectral en- 
ergy distribution (SED) built from their photometric datcQ 
to SED tem plates obtained from stella r population mod 
(such 



els 



Bruzual fc Chariot 



- - Mi3; ... 

(lLai et all feOQSl; IGawiser et al.l 120071: iNilsson et al.l 12007 . 



2003 ), hereafter BC03) 



| 2009l;lYabe et al Il200d; iFinkelstein et al.ll2009al . l2008l , l2009cl 
boiOal ; iMagdis et al.ll20ld ) . This procedure, in principle, en- 
ables the determination of age, dust attenuation, star forma- 
tion rate (and histories), metallicity, and stellar mass. This 
is because the SED obtained in stellar population models 
depends (among other) on all those parameters. However, 
in practice, this procedure has several limitations. For ex- 
ample, metallicity does not have a strong influence on the 
shape of the templates and, therefore, its determination from 
SED-fitting suffers from large uncertainties. On the other 
hand, the degeneracy between dust attenuation and age or 
age and star-formation history (SFH) produces that both 
parameters are difficult to constrain simultaneously. With 
a good wavelength coverage of the UV continuum and the 
Balmer break it is feasible to improve the determination 
of the SED-derived dust attenuation and stellar age. How- 
ever, the best determination of dust attenuation would still 
suffer from uncertainties and the only way to obtain accu- 
rate values is by employing direct FIR detections around 
the dust emission peak. Despite these caveats, many previ- 
ous works have analyzed the physical properties of LBGs at 
different redshifts by employing a SED fitting method since 
it is the only way to analyze their properties with large sam- 
ples of galaxies. At z ~ 5, LBGs have been reported to be 
much younger (<100Myr) and to have lower stellar masses 
(1O 9 M0) than their analogs at z ~ 2.0 — 3.0 in a similar rest- 
fram e UV luminosity range ijVerma et al.ll2007l : lYabe et aJj 
|2009| ; IHaberzettl et al.ll2012i ) . 

Most previous works focus their studies on LBGs which 
are located at z > 2 — 3, where the Lyman Break is shifted 
to the optical and can be sampled with filters in ground- 
based telescopes. At z<2, the Lyman Break is located in the 
UV and LBGs can only be found with observations from 
space , for example with GA LEX l|Burgarella et al" 1 120061 . 
l2007l ; IHaberzettl et al 1 12012| ). HST and UVIS filter chan- 
nels (|Hathi et al. 20121 ). or Swift satellite (|Basu-Zvch et al.l 
2011). In this work, we aim at analyzing the physical prop- 
erties of a sample of 1225 GALEX-selected LBGs at z ~ 1 
located in the COSMOS field by using data coming from 
the Advanced Large, Homogeneous Area M edium Band 
Reds hift Astronomical (ALHAMBRA) survey jMoles et all 
|2008| ). which covers the optical range with 20 medium-band 
filters (width about 300A) and the near-IR with the tra- 
ditional JHKs filters. The combination of the ALHAMBRA 
survey with observations in other wavelengths (GALEX and 



In this work we apply the term SED to refer to a set of pho- 
tometric points. However, it should be noted that SED is also 
applied to spectroscopic data in many works. Some of the limi- 
tations quoted here in the SED-fitting technique just applies to 
photometric-SEDs but not to spectroscopic data 



Lyman break and UV-selected galaxies at z ~ 1 3 



IRAC) allows an unprecedented coverage of the UV contin- 
uum and optical Balmer break which can disrupt some of 
the degeneracies outlined before and provide accurate re- 
sults for the physical properties obtained with SED-fitting 
methods. Since we study LBGs at z ~ 1, their observed 
fluxes are high enough so that the photometry has a good 
signal-to-noise ratio. These two facts (exceptional coverage 
of the SED and the high signal-to-noise ratio) are not usu- 
ally achieved at higher redshifts, which emphasizes the im- 
portance of studying intermediate-redshift LBGs utilizing 
the photometric coverage of the ALHAMBRA survey. 

The paper is organized as follows: In Section [2] we ex- 
plain the data sets employed both in the UV with GALEX 
observations and in the optical-to-near-IR with the AL- 
HAMBRA survey. In Section [3] we combine these two main 
data sets and build a general sample of UV-selected galax- 
ies at < 2 < 2. In that Section we also explain how we 
carry out the SED fits using BC03 templates for these UV- 
selected galaxies with the aim of obtaining their photomet- 
ric redshift, rest-frame UV luminosity, and other physical 
properties such as dust attenuation, age, and stellar mass. 
From this GALEX-selected galaxies we look among them 
the GALEX-selected LBGs and discuss their SED-derived 
physical properties in Section [5] The morphology and phys- 
ical sizes of the LBGs studied are analyzed in Section[6] and 
in Section [7] we show their location in the color-magnitude 
(CMD) diagram. In Section [8] we compare the properties 
of our GALEX-selected LBGs to those reported in previous 
works for LBGs at higher redshifts. Finally, in Section[9]the 
conclusions of this work are presented. 

Throughout this paper we assume a flat universe with 
(fi m ,fiA,/io) = (0. 3,0.7,0.7), and all magnitudes are listed 
in the AB system (|Oke fc Gunnlll983l '). 



2 DATA SETS 

On the UV side we use data coming from o bservations of 
the C OSMOS field with the GALEX satellite ijMartin et all 
120051 ) in both NUV and FUV bands as part of the Deep 
Imaging Survey (PI: D. Schiminovich) . GALEX catalogs 
were created by using the EM-algorithm ijGuillaume et all 
2006). This approach employs optical feel resolved cata- 
logs to resolve blended objects in the far and near UV. 
In the case of the COSMOS field, the prior optical photo- 
metric informations corresponds to a u*-band mosaic (and 
its SExtractor-derived catalog) based on CFHT-u* observa- 
tions of the COSMOS field. 

On the optical and near-IR side we use the ALHAM- 
BRA survey (|Moles et alj|2005l . I2OO8I ). which employs a set 
of 20 equal-width (~ 300 A) medium-band filters covering 
the optical range from 3500 to 9700Aplus the traditional 
JHKs broad-bands near-IR filters to observe a region of 4 
square degrees distributed into 8 different fields. The obser- 
vations have been obtained with the 3.5m telescope of the 
Calar Alto observatory using the wide-field camera LAICA 
in the optical and the OMEGA-2000 camera in the near-IR. 
Since we focus our work on COSMOS, we only consider the 
ALHAMBRA observations in that extragalactic field. The 
optical filters system adopted in ALHAMBRA was set with 
the aim of optimize t he output of the s urvey in terms of 
the photo-z accuracy (|Bem'tez et alJ koOQ). The simulations 



performed in iBem'tez et all (2009) relating the image depth, 
photo-z accuracy, and number of filters indicate that the fil- 
ter system of ALHAMBRA enable to get a photo-z precision, 
for normal SF galaxies, that is three times better than that 
for traditional 4-5 broadband filter sets. Additionally, the 
complementary usage of near-IR data improves the determi- 
nation of photometric redshifts. In this work we have used 
the catalogs coming from the Internal Data Release 3. The 
data reduction and the construction of the catalogs were car- 
ried out by the ALHAMBRA team. Since the ALHAMBRA 
survey performs observations in 23 filters, one must work 
with care when defining the detection of the objects. With 
the aim of not biasing the detection to any special kind of 
object, a special technique was employed. It is based on cre- 
ating a deep detection image built as the sum of a set of indi- 
vidual frames. These frames are those observed in the optical 
filters with the highest efficiencies and taken under the best 
observational conditions. This includes the filters centered 
between 457 and 829 nm. Then, the object detection is c ar- 
ried out by running SExtractor l|Bertin fc Arnoutsll 19961 ) in 
this auxiliary deep image. With the position of the detected 
objects in hand, their ALHAMBRA multi-wavelength pho- 
tometry is then obtained by running SExtractor in its stan- 
dard double image method applied to each of the filters. 
As result, the average depth (for 3<r detections) is 24.5 and 
22 magnitudes in the optical and near-IR, respectively. The 
characterization in the optical ran ge of the ALHAMBRA 
photo metric system can be found in Aparicio Villegas et alj 
l|2010l) and the near-IR number counts of one of the fields 
is presented in ICristobal-Hornillos et all (|2009r ). Details on 
the quality of the data, the reduction process, the depth, etc. 
will be published in Husillos et al. (in prep). We note that 
the characterization of ALHAMBRA filters (complete wave- 
length coverage with almost no overlapping filters) provides 
a SED that can be considered as a low resolution optical 
spectra of the observed sources. In this work, and owing to 
the wealth of data in many other wavelengths, we only ex- 
ploit ALHAMBRA data in the COSMOS field. It should be 
noted that the selection criteria of the LBGs that will be 
studied in this work are purely based on UV GALEX data, 
whereas the ALHAMBRA survey is used in the analysis of 
their SED-derived physical properties. 



3 GALEX AND ALHAMBRA DATA: 
UV-SELECTED GALAXIES AND SED 
FITTING 

The LBGs that will be studied in this work are taken from a 
catalog of UV-selected SF galaxies that is built by combin- 
ing the GALEX observations on the COSMOS field with the 
data coming from the ALHAMBRA survey. In order to com- 
bine GALEX and ALHAMBRA data we look for GALEX 
detections around 2" of the optical position of the sources 
in the ALHAMBRA survey and retain only those galaxies 
which have, at least, a detection in the NUV band. This 
produces a sample of UV-selected galaxies formed by 39734 
sources for which we have photometric information from the 
NUV to the near-IR. Since so far we have only imposed the 
detection in the NUV and optical-to-near-IR bands, it is ex- 
pected that these 39734 UV-selected galaxies have a wide 
range of both redshift and rest-frame UV luminosity. 
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3.1 SED fitting 

With the aim of both obtaining photometric redshifts and 
analyzing the physical properties of our sample of 39734 
UV-selected galaxies we use a SED fitting procedure to 
fit their observed fluxes with BC03 templates. In this pro- 
cess, the photometric SEDs are built by employing AUTO 
magnitudes. We build a set of BC03 templates associ- 
ated to different physical properties of galaxies by using 
the so ftware GALAXEV. In this process we adopt a ISalpeterl 
<ll955h initial mass function (IMF) distributing stars from 
0.1 to 100 Mq and select a constant value for metallicity 
of Z=0 .4Zp). Dust attenuati on is included in the templates 
via the lCalzetti et all (|2000T ) law and parametrized through 
the color excess in the stellar continuum, E S (B — V), for 
which we select values ranging from to 0.7 in steps of 0.05. 
We also include intergalac tic m edium absorption adopting 
the prescription of lMadaul (|l995h . Regarding SFR, we adopt 
time-constant models. In this case, different values of the 
SFR does not change the shape of the templates, and the 
SFR can be obtained by using the lKennicuttl (| 1998T ) calibra- 
tion: 



SFR 



UV, uncorrected 



[Moyr" 1 ] = 1.4 • lO^Lis 



(1) 



where L1500 is the rest-frame UV luminosity in 1500A, which 
is obtained from the best fitted template of each object by 
convolving it with a top-hat filter centered in 1500 A. It 
should be noted that, throughout the work, we distinguish 
between Luv defined in a v\j v way and L1500 considered in 
L„ units. The SFR derived from Equation [1] is uncorrected 
for the attenuation that dust produces in the rest-frame UV 
continuum. 

It should be noted that in this work we employ BC03 
templates associated to a time- independent SFR. Other kind 
of SFHs can be used, such as exponentially declining or com- 
posed by different burst of star formation. In the first case, 
the SFR is characterized by the decaying time-scale, tsfr, 
which would be another parameter to obtain in the SED fit- 
ting, increasing the degrees of freedom in the process. Distin- 
guishing between different kind of SFHs is very challenging 
even with a good photometric coverage of the SED of galax- 
ies. Therefore, the results reported in this work should be 
understood as those derived with that choice of the SFR, but 
different values of the SED-derived parameters might be ob- 
tained if other temporal variations of the SFH are assumed. 
The analysis of the differences in the SED-fitting results de- 
pending on the assumption of the SFH will be studied in 
Section [5T2l 

Once the dust attenuation is known for each source, we 
can obtain an estimation of the dust-corrected total SFR by 
correcting the observed L1500 from dust absorption. To do 
that, we take the observed L1500 and multiply it by the factor 
10°' 4Al5O °, where A1500 is the dust attenuation in 1500A, 
which can b e obtained from the SED-derived E a (B — V) 
assuming the lCalzetti et al. (2000) law. Once both age and 
dust-corrected total SFR are known for each source, and 
according to the assumed time- independent SFH, the stellar 
mass can be obtained from the product of both quantities. 

In this work we also aim at analyzing the UV contin- 
uum slope, ft, of our UV-selected galaxies. This parameter 
is important in the study of how galaxies built up since 



it is related to age, metallicity, stellar IMF, and most im- 
portantly, dust attenuation. Furthermore, UV colors seem 
to present a relation with UV lum inosities of SF galaxies 
jBouwens et al.ll2009l . l2010al . | 2011al ) and are relatively easy 
to measure than optical rest-frame colors in high-redshift 
galaxies, for which IRAC detections are mandatory. At our 
redshift range, the combination of the GALEX photometry 
and the bluest optical bands of the ALH AMBRA survey pro- 
vides a good sampling of the UV continuum in a wide range 
of redshifts, as it can be seen in the SED-fitting results shown 
in Figure [T] Different works employ different methods to 
obtain the UV continuum slope of galaxies at different red- 
shifts ranges, being the most popular and traditionally used 
that in where /3 is quantified by using two broadband fil- 
ters which sample t wo zones of the ob s erved UV cont i nuum 
jKong et all |2004|; lHathi et al.l |2008|; iMeurer et all Il997l; 
Overzier et al.l2008l;|Bouwens et al.l2010bl ; lFinkelstein et al.l 
2010bl ; lDunlop et alj|2012h . In other works /3 is obtained by 



using a power-law fit to the observed fluxes of the studied 
galaxies, using filters t hat sample the same z one of the SED 
at different redshifts (|Bouwens et all 1201 lbl ). In our work, 
we obtain /3 for each galaxy directly from its best-fit tem- 
plate by fitting a power law in the form fx ~ A' 3 to their 
UV continuum. In this process we employ the rest-frame 
wavelength range between 1800A < A < 3000A. This ap- 
proach has the advantage of using all the available fluxes of 
each source, resulting in more robust S/N determinations, 
whereas the results based on pairs of colors can suffer from 
large uncertainties due to photometric out liers. The method 
emplo yed here is similar to that used by iFinkelstein et ail 
(|20 1 IT ) in their study of the redshift evolution of the UV 
continuum slope from z ~ 8 to z ~ 4. In that work, they 
present some illustrative examples showing the differences 
in the UV continuum slope when using different techniques. 

The fits are carried out by using the Z urich Extragalac- 
tic B ayesian Redshift Analyzer (ZEBRA, Fe ldmann et al.l 
2006) code which, in its maximum-likelihood mode, employs 
a x 2 minimization algorithm over the templates to find the 
one which best fits the observed SED of each input object. 
Each fitted template has a characteristic value of the x 2 . 
Here we define the reduced x 2 1 Xr, as the ratio between 
the x 2 associated to each fitted template and the number 
of filters minus one employed in the SED-fitting procedure 
for each source, Xr — X 2 /(N — 1). From a visual inspec- 
tion of the SED-fitting results, we consider that the BC03 
templates truly represent the observed SED for each galaxy 
when Xr < 10 (see some examples of x 2 values and the 
quality of the fittings in Figure [TJ. Imposing x 2 < 10 to the 
fittings of the whole sample of 39754 UV-selected galaxies, 
we end up with a robust sample of 35810 galaxies. From now 
on, only sources with Xr < 10 are considered. 

In a SED-fitting procedure, the reliability of the results 
i.e. the similarity between the observed SED and that rep- 
resented by its best fitted template, is related to the x 2 
value of the fits. As indicated above, a visual inspection of 
the SED-fitting results of our UV-selected galaxies makes 
us to conclude that the SED fittings, and therefore their 
results, are reliable when Xr < 10, which is the threshold 
value that we adopt here. It should be considered that the 
Xr values depend upon the observed fluxes and also upon 
their photometric uncertainties. In this way, if a galaxy has 
a photometry with high photometric errors, the Xr might be 
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low even when its best fitted template does not represent its 
observed SED properly. A low value of the Xr can be due to 
either a good SED-fitting or to a 'not-so-good' SED-fitting 
with a photometry with relatively high uncertainties. Thus, 
we should check the typical photometric errors of the AL- 
HAMBRA photometry of our sources to analyze if the low 
values of Xr are due to a truly accurate fit or if they are the 
consequence of high photometric errors. As an example, we 
show in the left plot of Figure [2] the photometric errors in 
the ALHAMBRA filter centered in 613 nm of our sample of 
35810 UV-selected galaxies as a function of their observed 
magnitude in the same band. As expected, the photometric 
errors increase with the observed magnitude. If we consider 
that a fit is reliable for galaxies with typical photometric 
errors below 0.4 mag, we can only trust those SED-fitting 
results for galaxies typically brighter than about 24 mag. In 
the right plot of Figure [2] we represent the observed magni- 
tudes of the GALEX-selected LBGs that will be studied in 
this work. It can be seen that most of them meet the pre- 
vious criterion and, therefore, we can consider that the low 
values of Xr are due to good SED fits rather than to high 
photometric errors. 
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3.2 Photometric redshifts and their accuracy 

The good coverage of the observed UV-to-near-IR SED of 
galaxies provided by the ALHAMBRA survey in combina- 
tion with GALEX data is expected to give accurate deter- 
mination of photometric redshifts (z p hot) at the expected 
redshift range of GALEX-selected LBGs, i.e. z ~ 1. This is 
due to the fact that at z ~ 1, GALEX+ ALHAMBRA data 
covers the rest-frame UV continuum and the Balmer break, 
which are two of the most important features to fit in the 
SED of galaxies to determine their photometric redshifts. 
To check the accuracy of the photometric redshifts, we com- 
pare in Figure [3] the obtained z p hot with spectroscopic red- 
shifts (zspec) for those sources in the whole sample of 35810 
UV-selected galaxies with ALHAMBRA measurements and 
Xr < 10 in the SED-fitting result s which have avai lable spec- 
tra from the zCOSMOS survey l|Lillv et alj|2007h . Here we 
define the accuracy of the z v hot a s ctaz = Az/( 1 + z B vec ), 
with Az = \z p hot ~ Zspec]' See also lMatute et al] l|2012h for 
a discussion of the photometric redshift accuracy of AL- 
HAMBRA survey. It can be seen in Figure [3] that within 
the redshift range 0.8 < z < 1.2 there is a good agree- 
ment between the photometric and spectroscopic redshifts, 
being ctaz less than 0.05 for most galaxies (see horizontal 
dashed lines). Considering a median redshift of 0.75, that 
value of (JAz translates into a maximum redshift uncertainty 
of Az ~ 0.04. It should be noted that this accuracy only 
applies to galaxies which are as bright as the sources in 
the spectroscopic survey. All the galaxies with spectroscopic 
redshift from zCOSMOS survey shown in Figure [3] have r- 
band observed magnitudes typically brighter than 23.5 mag 
and, therefore, the reliability of the photometric redshifts 
can be guaranteed up to that limit. From now on, photo- 
metric redshifts are used for the sources without zCOSMOS 
counterpart. For those sources with a zCOSMOS spectrum, 
we use the spectroscopic redshifts. 



Figure 5. Location of our GALEX-selected LBGs in a color- 
color diagram. The window enclosed by the orange lines is the 
selection region for LBGs at 0.6 < z < 1.4 with NUV < 23.75. 
The subsample of our GALEX-selected LBGs which satisfies 
NUV < 23.75 are represented with red filled dots, while the re- 
maining fainter LGBs are plotted with black symbols. 



4 UV-SELECTED GALAXIES ATZ~1 

Once we have discussed the selection and built of the whole 
sample of 35810 UV-selected galaxies with good SED-fitting 
and have obtained their photometric redshifts, we are in 
position to explain the selection of our LBGs at z ~ 1. We 
sample the Lyman break of galaxies at z ~ 1 by employing 
their emission in the FUV and NUV GALEX filters. 

As it was commented in Section [TJ the choice of the 
red-ward and blue-ward filters determines the wavelength 
where the Lyman break is located and, therefore, the range 
of redshifts of the selected galaxies. In this work we aim 
at analyzing the physical properties of those LBGs whose 
Lyman break is located between the GALEX FUV and NUV 
filters, which are centered at 1528A and 227lA (in terms of 
their effective wavelengths) and have bandwidths of 1344- 
1786 A and 1771-2831 A, respectively. These wavelengths 
imply that the redshifts of these GALEX-selected LBGs are 
expected to be around z ~ 0.95 considering an intermediate 
wavelength of 1780 A between the two filters. 

In order to formulate an analytic selection criteria to 
segregate our LBGs and that could be applied to the ob- 
served magnitudes of the sample of UV-selected galaxies 
previously built, we use a large set of BC03 templates associ- 
ated to a metallicity of Z = OAZq, constant star formation 
rate, and different values of age and dust attenuation. We 
study the location of those templates in a color-color-redshift 
space. To do that, each template is redshifted from z — 
up to 2 = 2, then we ap ply the correspondi ng absorption 
in the IGM following the lMadau et all l|l996l ) prescription, 
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Figure 1. Examples of SED-fitting results with GALEX+ALHAMBRA data for nine GALEX-selected LBGs randomly selected within 
the whole sample. Blue points are the observed GALEX and ALHAMBRA fluxes and the black curves are the BC03 templates that fit 
the photometry of each object best. In each plot we indicate the SED-derived redshift, age, dust attenuation, and reduced \ 2 °f t ne 
fittings. It can be seen that the combination of GALEX and ALHAMBRA data provides a very good sampling of the rest-frame UV 
continuum and Balmer break of our UV-selected galaxies. 



and obtain the FUV — NUV synthetic observed colors by 
convolving the templates with the transmission curves of the 
GALEX filters. The result is shown in Figure [4] As a general 
trend and as it could be expected by the location of the Ly- 
man break as a function of redshift, the FUV — NUV color 
increases with redshift up to z ~ 1. Looking at the different 
tracks represented in Figure [4J we decide to impose a color 
cut of 1.5 (red dashed horizontal line) and, therefore, our 
color selection for GALEX-selected LBGs is: 

FUV - NUV > 1.5 (2) 
It is important to note that the application of this crite- 



rion requires the detection of each galaxy in both FUV and 
NUV channels. Left plot of Figure [4] indicates that imposing 
such a color selection criterion we segregate galaxies located 
at z > 0.7 (this threshold is represented by the red dashed 
vertical line). However, at 0.7 ^ z ^ 0.8, the FUV flux is 
importantly affected by the photons of the Lyman contin- 
uum (Lyc), which are those UV photons whose wavelengths 
are higher that the wavelength of the Lyman break at each 
redshift. Therefore, if we really want to sample the Lyman 
break between the FUV and NUV without significant con- 
tamination of the Lyc in the FUV filter we have to limit the 
redshift of the galaxies at z ^ 0.8. This situation is schema- 
tized in the right plot of Figure 2] It can be seen that the 
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Figure 2. Left; Photometric errors against the observed magnitude in the ALHAMBRA filter centered in 613 nm for our sample of 
UV-selected galaxies with ALHAMBRA measurements and Xr < 10 m the SED-fitting results. Red curve corresponds to the median 
value of the error distribution for each value of the observed magnitude. Horizontal lines represent photometric errors of 0.1 and 0.2 
mag. Right: Distribution of the observed magnitude in the ALHAMBRA filter centered in 613 nm for LBGs and a sample of UV-selected 
galaxies with ALHAMBRA counterparts at the same redshift range than LBGs. Histograms have been normalized to their maxima in 
order to clarify the representation. 
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Figure 3. Accuracy of the photometric redshift determination with the combination of GALEX and ALHAMBRA data for the whole 
sample of 35810 UV-selected galaxies with ALHAMBRA me asurements and % j < 10 m the SED-fitting results. In these plots, galaxies 
with available spectroscopic redshift from zCOSMOS survey dLillv et al.ll2007Tl are considered. Red dots represent our GALEX-selected 
LBGs and grey open triangles are the remaining galaxies in the sample. In the left plot, vertical and horizontal dashed straight lines 
represent the photometric redshift locus where most GALEX-selected LBGs are expected to be located according to their UV color 
selection. In the right plot, vertical dashed straight lines represent the photometric redshift locus where most GALEX-selected LBGs are 
expected to be located according to their UV color selection. The horizontal dashed straight lines indicates the values of the photometric 
redshift accuracy, = Az/(1 + z spec ), equal to ±0.05. The horizontal continuous straight line represents where spectroscopic and 

photometric redshifts would agree. 
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Figure 4. Color selection of the GALEX-selected LBGs studied in this work. Left plot shows the synthetic FUV-NUV color tracks as 
a function of redshift according to a set of BC03 stellar population templates associated to different values of age and dust attenuation. 
The horizontal dashed red line indicates the color cut employed in this work, which according to the tracks shown in black is expected 
to select galaxies at z J 0.7 (vertical red dashed line). In the left plot we represent with orange and red curves the transmission of the 
FUV and NUV GALEX filters, respectively. We also show the location of the Lyman break at z = 0.8 with a vertical dashed black line 
and the BC03 stellar population templates associated to three different values of age and dust attenuation (black curves). It can be seen 
that, although the color selection shown in the left plot is expected to segregate galaxies at z ^ 0.7, the Lyman break has not almost 
completely passed the FUV channel until z ~ 0.8. 



Lyman break at z = 0.8 is located almost out of the FUV 
filter and therefore the contamination in that filter from Lyc 
photons is negligible. At lower redshifts, there would be a 
significant contribution of the Lyc photons in the FUV filter 
and, therefore, we would not be truly sampling the Lyman 
break, but the Lyman continuum. It should be noted that, 
as the dropout technique is a method to find high-redshift 
galaxies, the sampling of the Lyc continuum with both FUV 
and NUV filters would be a perfectly valid technique to look 
for galaxies at 0.6 ^ z ^ 0.8. However, we decided to select 
only galaxies in which the Lyman break is located between 
the FUV and NUV filters with low contamination of the Lyc 
to the FUV channel and, the refore, consider galaxies only 
at z ^ 0.8. As a comparison, iBurgarella et all (120071) con- 



sider only galaxies with 0.9 ^ z ^ 1.4, whereas Mosleh et al.l 
|201ll ) retain all the gal axies with 0.6 ^ z ^ 1.4 following 
the selection criterion of lBarger" et al.1 l|2008l ). 

Then, we define, as a first approximation, GALEX- 
selected LBGs as those galaxies which are detected in both 
FUV and NUV channels and whose fluxes in each band meet 
the Equation [2] and are located at z ^ 0.8. This sample is 
formed by 475 galaxies. It is worth remarking that there is a 
difference between the selection criterion that we apply here, 
and those applied to look for high-redshift LBGs. At z > 2, 
LBGs are usually found employing not only the difference 
in color which characterizes the Lyman Break (Equation [2] 
in our case) but also a d ifference in color in redder wave- 
lengt hs (see for example ISteidel et al.l 120031 ; iMadau et al.l 
1 19961 ) for ruling out interlopers, i.e. galaxies at other red- 



shifts that have a large difference in color between two con- 
tiguous broad-bands filters but not due to the Lyman break. 
This has to the done because at high redshift the photomet- 
ric redshifts tend to suffer from large uncertainties and then, 
it is not possible to select galaxies at a specific redshift range 
using z-phot. However, in our case, as it was commented in 
Section \3. 21 and shown in Figure [3l we have accurate values 
of the photometric redshift for our UV-selected galaxies at 
z ~ 1 a nd, therefore , that supplementary condition is not 
needed. IBarger et all (|200St ) select LBGs at 0.6 s; z 1.4 
by employing a double color selection criterion combining 
FUV - NUV and NUV - U. If we limit our sample in NUV 
magnitude to their same limit, NUV < 23.75, all but one 
of our GALEX- selected LBGs satisf y the double color selec- 
tion criterion of lBarger et al.l (|2008T l. This is schematized in 
Figure [5] The U-band data for the galaxies in the plot have 
been taken from the broad-band ph otometric catalog in the 
COSMOS field (|Capak et al.ll2007l ). Conversely, if we trust 
the photometric redshifts obtained from the combination of 
GALEX and AL HAMBRA da t a we find that with the selec- 
tion criterion of lBarger et al.l (|2008l ) we would miss a pop- 
ulation of GALEX-selected LBGs fainter that NUV = 23.5 
mag whose NUV — U color are typically redder that those 
GALEX-selected LBGs with NUV < 23.5 mag. 

The application of Equation [2] requires the detection 
of the galaxies in both the FUV and NUV GALEX chan- 
nels so that the amplitude of the break can be measured. 
However, it can also occur that a galaxy has such a strong 
Lyman break that is is detected in the NUV but its FUV 
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Figure 6. Distribution of the FUV — NUV color for Xray- 
detected galaxies at 0.8 < z < 1.2 with GALEX and ALHAM- 
BRA counterparts (black histogram). Red histogram represents 
the distribution of the FUV — NUV color of the galaxies at 
0.8 z ^1.2 spectroscopically classified as AGNs in lCowie et alj 
via emission line diagnosis. The vertical dashed line in- 
dicates the color threshold for selecting LBGs in this work. In 
this plot, only galaxies with detection in both FUV and NUV 
channels are included. 



magnitude is fainter than the FUV limiting magnitude and, 
therefore, is not detected in the FUV channel. In these cases, 
although Equation [2] cannot be strictly applied, these galax- 
ies are also LBGs due to their strong Lyman break between 
the FUV and NUV filters. Following this idea, we also se- 
lect those galaxies with are detected in the NUV channel 
but are undetected in the FUV one. In order to include 
these FUV-undetected galaxies into the sample of GALEX- 
selected LBGs (i.e. to be confident that they truly have a 
strong Lyman break) we have to confirm that they do have 
a strong Lyman break which causes the non-detection in the 
FUV channel. FUV observations in the COSMOS field have 
a limiting magnitude of FUV ~ 27 mag. Galaxies brighter 
than that value in the wavelength range covered by the FUV 
filters should have been detected. Since we select galaxies 
with a FUV — NUV color cut of 1.5, that limiting magni- 
tude would imply a limit of 25.5 mag in the NUV channel. 
To be conservative, we choose a limit in the NUV channel 
of 25 mag and include in the previous sample of LBGs those 
galaxies that are brighter than 25 mag in the NUV channel, 
are at 0.8 z 1.2, and are undetected in the FUV chan- 
nel. With this sample included, we end up with a sample 
of 1246 GALEX-selected LBGs. A visual inspection of the 
galaxies with available ACS information (see Section [5] for 
more details) reveals that the contamination due to the low 
spatial resolution of the GALEX images is lower than 5%. 



4.1 X-ray counterparts and AGN contamination 

The Lyman break selection technique segregates galaxies 
which are bright enough in the UV continuum so that the 
Lyman break can be sampled by using two broad band fil- 
ters. In principle, the UV emission detected in LBGs can 
come from hot stars or AGN accretion disks. As we are only 
interested in those LBGs which are SF galaxies, we need to 
rule out the AGN contribution. This can be done either with 
rest-frame UV/optical spectroscopy or by means of their 
X-ray emission. Since we do not have spectroscopy for all 
the galaxies selected in the previous section, we have to ex- 
plore the possible AGN contamination using X-ray observa- 
tions. To this aim , we look for CHANDRA X-ray detecti ons 
jElvis et all 120091 ) around 3" (|Povic et all 120091 , l2012h of 
the ALHAMBRA-based spatial coordinates of our GALEX- 
selected LBGs. The area where our LBGs are located is al- 
most totally covered by the CHANDRA footprint. We find 
that only 21 LBGs are detected in X-rays and, therefore, 
likely have an AGN na ture. By using the catalog of AGNs 
in the COSMOS field l|Salvato et al.ll201ll ) we do not find 
any extra AGN identification. The only AGN-LBGs repre- 
sents an AGN contamination of about 2%. 

Figure [6] represents (black histogram) the distribution 
of FUV — NUV colors of all the galaxies with detection 
in FUV and NUV channels which are detected in X-rays, 
have GALEX and ALHAMBRA counterparts, and are at 
0.8 ^ z ^ 1.2 (the redshift range where our GALEX-selected 
LBGs are located). It can be seen that most X-ray-detected 
galaxies and, therefore, galaxies with an AGN nature, have 
FUV — NUV colors below the color threshold utilized in this 
work for selecting LBGs (see Equation]!?]). This UV color dis- 
tribution for AGNs at z ~ 1 explains the low percentage of 
AGNs among the GALEX-selected LBGs. We also plot in 
Figure [S] the FUV — NUV color distribution of the galaxies 
spectros copically clas s ified as AGNs via emission line diag- 
nosis in ICowie et al,l (|20101 ) and which are located at the 
same redshift range of our GALEX-selected LBGs. It can 
be seen again that most AGNs at z ~ 1 have FUV — NUV 
colors below the color threshold considered in this work for 
selecting LBGs, reinforcing the fact that the LBG color se- 
lection does not tend to segregate galaxies with an AGN 
nature. 

Low values of the AGN contribution in sa mples of LBGs 
at diff erent redshifts have been also reported. iLehmer et al.l 
(|2005h found AGN fractions of 1.2%, 0.4%, 0.3% and 0.4% 
in their sa mple of U-, B435 - , Vjiog -, and i775-dropouts, re- 
spectively. iBasu-Zvch et all |2011) reported an AGN frac- 
tion for their sample of LBGs at 0.5 < z < 2.0 of 5%-6%, and 
iNandra et all (2002) found an AGN contribution of about 
3% in their sample of LBGs at z ~ 3. In the subsequent anal- 
ysis we do not take into consideration the GALEX-selected 
LBGs with active nuclei. Therefore, we end up with a sam- 
ple of 1225 SF GALEX-selected LBGs. This is the largest 
sample of LBGs studied at z ~ 1 so far. Figure [1] shows 
the UV-to-near-IR SED of nine GALEX-selected LBGs in 
our final sample. This small subsample is representative of 
the whole sample of LBGs. It can be clearly seen that the 
combination of GALEX and ALHAMBRA provides an ex- 
cellent coverage of the rest-frame UV continuum, Balmer 
break, and near-IR SEDs of these galaxies. 

Figure [7] represents the distributions of the apparent 
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Figure 7. Distribution of the NUV (blue histogram), ALH-706 
(red histogram), and Ks (red histogram) apparent magnitudes for 
our GALEX-selected LBGs at z ~ 1. 

brightness of our GALEX-selected LBGs in the NUV chan- 
nel, optical ALH-706 ALHAMBRA filter, and in the near-IR 
Ks band. These distributions should be taken into account 
for comparing our results for our LBGs with those reported 
in other works that employ different photometric informa- 
tion. 

4.2 High-redshift analogues 

In this work, we also aim at comparing the SED-derived 
physical properties of LBGs at different redshifts. At high- 
redshift, i.e. z > 3, and due to the observational bias pro- 
duced by the use of magnitude-limited catalogs, LBGs at 
high redshift tend to be intrinsically more luminous than 
those selected in the present work. Therefore, if we want to 
compare LBGs at different redshifts and, therefore, compare 
galaxies which are selected following similar selection crite- 
ria we must limit our sample to the same rest-frame UV 
luminosity range than that for high-redshift LBGs, which 
is typically logLuv/i© *5 10.2. In this sense, we define 
UV-bright LBGs as those LBGs at 0.8 < z < 1.2 which 
have logLuv/L© 10.2. This sub-sample is formed by 181 
galaxies. 

4.3 UV-faint galaxies 

The Lyman break technique segregates, at each redshift, 
galaxies which tend to be UV-bright. In our case, most 
GALEX-selected LBGs at 0.8 < z < 1.2 have NUV mag- 
nitudes brighter than 25. In that redshift range, there are 
many UV-detected SF galaxies which are intrinsically fainter 
than the GALEX-selected LBGs and, as a consequence, the 
could not been selected using the dropout criterion. All these 
galaxies which tend to be intrinsically UV fainter than the 



GALEX-selected LBGs and to which the dropout criterion 
could not be applied, will be called UV-faint galaxies. This 
sample will not be studied in the present work since we aim 
here to study the physical properties of GALEX-selected 
LBGs. However, it will be used in an incoming work (Oteo 
et al. in prep.) where FIR observations will be used to con- 
strain the FIR SED of both GALEX-selected LBGs and 
UV-faint galaxies. In that case, the comparison between 
FIR-detected LBGs and UV-faint galaxies will help to un- 
derstand the galaxies which are selected under the dropout 
technique in opposition to other UV-fainter SF galaxies and 
to place LBGs in a more general scenario of SF galaxies at 
0-8 < z < 1.2. 



5 SED-DERIVED STELLAR POPULATIONS 

5.1 Physical properties of LBGs at z ~ 1 

Figure [8] shows with red shaded histograms the distributions 
of photometric/spectroscopic redshift, rest-frame UV lumi- 
nosities, age, dust attenuation, dust-corrected total SFR, 
stellar mass, and UV continuum slope for our GALEX- 
selected LBGs. It can be seen that, as a consequence of 
their selection criterion, they are located at 0.8 < z < 1.2 
and have rest-frame UV luminosities log (Luv/i©) > 9.6. 
It can be seen that GALEX-selected LBGs at 0.8 < z < 1.2 
are young objects, with ages mainly distributed within < 
Age[Myr] < 500 with a median value of 341 Myr. Their dust 
attenuation parametrized by the color excess in the stellar 
continuum E S (B — V) varies mainly from 0.0 to 0.4, with 
a median value of 0.20. This indicates that LBGs at z ~ 1 
are galaxies with moderate dust attenuation. Since they are 
UV-bright galaxies, their dust-uncorrected SFR are also rel- 
atively high, ranging typically between 1 and 6 M©yr _1 , 
with a median value of 2.OM0yr _1 . Their stellar masses 
span within 9 < log(Af*/M0) < 11, with a median value 
of log (M*/Mq) — 9.74. LBGs are also blue galaxies, with a 
UV continuum whose slope varies mostly between -2.5 and 0, 
with a median value of/? = —1.53. Finally, GALEX-selected 
LBGs have high values of the dust-corrected total SFR with 
a median value of 46.4 Moyr -1 . 

In its maximum likelihood mode, ZEBRA gives not only 
the best-fitted templates but also the probability that any of 
the others non-best-fitted templates can represent the pho- 
tometric SED of a given galaxy. This probability can be 
used for deriving the uncertainties of the SED-derived pa- 
rameters. To this aim, we define the weighted average (IfA) 
value of each parameter as: W A = Pifi/N, where Pi is 
the probability that a given template, i, can represent the 
observed SED of a given galaxy, fi is the value of one of 
the physical properties associated to that template, and N 
is the number of templates. For a given galaxy, if the best 
fitted template has much higher probability to represents its 
observed SED than any of the others templates, the average 
weighted value would be quite similar to the best fitted one 
and the uncertainty in the parameter should be low. On the 
other hand, if several templates associated to very differ- 
ent values of a given parameter have similar probability to 
represent the observed SED of a given galaxy, the weighted 
value would be dissimilar to the best-fitted one and the un- 
certainty of that parameter should be high. Following this 
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Figure 8. From left to right and from top to bottom, distributions of redshift (photometric or spectroscopic), rest-frame UV luminosity, 
and SED-derived age, dust attenuation, dust-corrected total SFR, stellar mass, and UV continuum slope for our GALEX-selected LBGs. 
BC03 templates associated to a constant SFR and metallicity Z = 0.4Zq are considered in the SED fits. Histograms have been normalized 
to their maxima. 



idea, we define the uncertainty of a given parameter as the 
difference between the best-fitted value and its correspond- 
ing weighted value. 

The physical parameters intrinsically related to the 
BC03 templates considered in this work for deriving the 
properties of the studied galaxies are the stellar age and the 
dust attenuation. The UV-derived dust-uncorrected SFR is 
obtained from the normalization of each observed template 
to the observed SED and the stellar mass is obtained from 
the already derived values of age, dust attenuation and dust- 
uncorrected SFR. This way, the procedure outlined above for 
deriving the uncertainties can be applied to age and dust at- 
tenuation. The median value of the uncertainties of age and 
dust attenuation for our LBGs are Aage = 390 Myr and 
AE S (B — V) — 0.1, respectively. The typical uncertainty of 
the SED-derived age is of the same order that the median 
age of our GALEX-selected LBGs. This implies that, even 
with the exceptional photometric coverage of the ALHAM- 
BRA survey, which samples quite well the Balmer break of 
SF galaxies at z ~ 1, the age is a parameter difficult to 
determine accurately with an SED-fitting procedure. 

The uncertainty in the SFR are directly related to the 
uncertainties of the rest-frame UV luminosity, related at the 
same time to the photometric uncertainties of the sources. 
The rest-frame 1500 A is sampled by the NUV filter at z ~ 1. 
The typical uncertainties of the NUV magnitudes for the 
galaxies in our sample is 0.15 This value translates to a 
luminosity uncertainty of A log (Ltjv/Lq) ~ 0.06 and an 
UV-derived SFR uncertainty of ASFR = 0.15Af Q /yr. 

INilsson et all l|201lf ) studied the optical SED-derived 
physical properties o f a sample of 15 GALEX -selected LBGs 
at z ~ 1, taken from lBurgarella et al.l (|2007l ). by fitting with 
BC03 templates their photometric points derived from ACS 
slitless grism spectra. They obtain similar distributions for 
dust attenuation and stellar mass than those reported in the 



present work, but th eir ages tend to be lower. In their work, 
INilsson et al.l l|201ll ) employed BC03 templates which are a 
combination of two SSP models instead of the models with 
constant SFR assumed in the present work and the ages re- 
ported are the youngest ones associated to the two SSP mod- 
els used in the fittings. This i s likel y the reason why the age 
distrib ution of Nilsson et al. (|201lh is shifted towards lower 
values. iBasu-Zvch et al.l (|201 W studied the physical proper- 
ties of a sample of 50 LBGs at 0.5 < z < 2 selected with the 
SWIFT ultraviolet/optical telescope. They reported lower 
values of stellar masses, (log iW* /M®) = 9.4±0.6, than those 
obtained here and with a similar dust attenuation distribu- 
tion. In this case, the differences are likely due to the fact 
that their galaxies spans within a redshift range which ex- 
pands towards higher values caused by the different filters 
used to segregate the galaxies. 

It is worth noting that the calculation of the dust- 
corrected total SFR has been obtained from a SED-fitting 
in the UV-to-near-IR regime by employing the SED-derived 
dust attenuation E a (B — V). However, this is not the best 
way of obtaining the dust attenuation, and consequently, the 
total SFR in galaxies. The most accurate procedure is uti- 
lizing direct measurements in the UV and in the FIR, since 
they probe the UV emission and the FIR emission of the 
light absorbed by dust in the UV, respectively. Depending 
on redshift, rest-frame UV luminosity, total IR luminosity, 
and kind of galaxy (LAE, LBG, SMG, ...), the SED-derived 
dust attenuation can be similar, overestimated, or underes- 
timated in comparison to t hat derived by emp loying direct 
UV and FIR measurements. IWuvts et al.l (|201ll ) analyze the 
relation between both quantities as a function of redshift by 
using direct measurements in the FIR from PACS photome- 
try. They use BC03 templates built by considering exponen- 
tial decreasing star formation histories and therefore their 
SED-derived dust attenuations might differ from ours, ob- 



12 I. Otto et al. 



CD 
M 



CD 

ao 
<^ 

BJO 

o 



0.25 0.53 0.B2 1.10 1.38 1.67 1.95 
3 

2 
1 




-1 
-2 
1.0 



• . : i »«:"*•?! 







1.5 2.0 2.5 3.0 3.5 4.0 
Log(Age 

constant 

[Myr]) 



<d 

< 



V 

<^ 

ao 
o 



0.245 0.531 0.816 1.102 1.388 1.673 1.959 

3 







1 
2 
1.0 



1.5 2.0 2.5 3.0 3.5 4.0 
Log(Age 

constant 

[Myr]) 



CD 
60 
< 



0) 
6fl 
<_ 

&0 
O 



0.40 0.64 0.87 1.11 

3 t 
2 
1 



1 r 

-2 L 



1.34 1.58 1.81 











cd 



CD 
M 



1.0 1.5 2.0 2.5 3.0 3.5 4.0ac 
Log(Age 

constant 

[Myr]) 3 




1.5 2.0 2.5 3.0 3.5 
Log(Age 

constant 

[Myr]) 



0.78 1.41 2.04 2.67 3.30 3.93 4.56 



1 

-2 

1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Log(Age 

constant 

[Myr]) 



0.90 1.02 1.14 1.26 1.39 1.51 1.63 



0.94 1.02 1.09 1.16 1.23 1.30 1.38 S 



0.9511.0171.0821.1471.2121.2771.342 



CD 
SO 
< 



CD 
M 
<^ 

^ao 
o 



3 
2 

1 



1 

-2 
1.0 





















1 1 '"' 







CD 
M 



CD 

ac 



1.5 2.0 2.5 3.0 3.5 4.0^ 



Log(Age c 



[Myr]) 



3 
2 

1 



-1 
-2 
1.0 



CD 

< 



CD 

< 



1.5 2.0 2.5 3.0 3.5 4.0g» 



Log(Age c 



[Myr]) 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Log(Age 

constant 

[Myr]) 



Figure 9. Differences between the SED-derived age when performing SED fittings with BC03 templates associated to constant SFR 
and BC03 templates associated to exponentially declining SFH with finite values of the SFH time scale. The color of the points in each 
figure is related to the ratio in the Xr values between the SED-fitting results when considering the different kinds of SFHs. The values 
of such Xr ratios corresponding to each color are indicated by the color bars. Values close to one indicate that the templates associated 
to different SFHs fit the photometry with the same accuracy and, therefore, it is not possible to distinguish between different kinds of 
SFHs. Each plot is associated to one value of the SFH time scale, as indicated in each vertical axis. 



tained with BC03 templates associated to constant SFR. 
Despite this, they find that at z ~ 1 and for modest levels 
of star formation both estimations agree well in the median 
values although with some scatter. This could be the case of 
IR-bright and UV-obscured sources. Actually, some LBGs at 
different redshifts have been reported to be IR-bright galax- 
ies. O ne of the most recent works is due to iBurgarella et"al] 
(2011), who report direct detections in the FIR by using 
SPIRE- Herschel data, probing directly the dust emission of 
LBGs at z ~ 1.0. The derived IR luminosities of those galax- 
ies place them into the LIRG/(U)LIRG regime. Some other 
works have also reported FIR/Sub-mm detection in LBGs, 
which also reveals the IR-brig ht nature of a subsample of 
them ijChapman fc Casevlliooa ). 



5.2 Dependance on the assumed SFH 

The physical properties of our GALEX-selected LBGs 
that we have discussed in the previous section were de- 
rived by carrying out a SED fitting procedure with BC03 
templates built by assuming a constant SFR. Different 
works assume different SFH for the galaxies under study 
llNilsson et al.l201ll;lBasu-Zvch et al.l201ll ; lHaberzettl et all 
120121 : lHathi et al.ll2012f ). One of the most used is that where 
the SFH varies exponentially with time. Analytically, it can 
be described by SFH oc exp (— i/rsFii), where tsfh is the 
SFH time scale. It should be noted that the constant SFR 
is a particular case of the exponentially declining SFH when 
the SFH time-scale tends to infinity. The usage of exponen- 
tially declining SFHs entails the fitting of the SFH time-scale 
in the SED-fitting procedure, which increase the number of 
degrees of freedom. Furthermore, there is a degeneracy be- 
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Figure 10. Differences between the SED-derived stellar mass when performing SED fittings with BC03 templates associated to constant 
SFR and BC03 templates associated to exponentially declining SFH with finite values of the SFH time scale. The color of the points 
in each figure are related to the ratio in the Xr values between the SED-fitting results when considering the different kinds of SFHs. 
The values of such Xr ratios corresponding to each color are indicated by the color bars. Values close to one indicate that the templates 
associated to different SFHs fit the photometry with the same accuracy and, therefore, it is not possible to distinguish between different 
kinds of SFHs. Each plot is associated to one value of the SFH time scale, as indicated in each vertical axis. 



tween the SED-derived age and the SFH time scale that 
produces that different combinations of age and tsfh give 
templates with the same shape and, therefore, non distin- 
guishable with photometric SEDs. In this Section we ana- 
lyze the differences in the SED-derived physical properties 
for our GALEX-selected LBGs when assuming SFH with 
different values of the SFH time-scale. Furthermore, we also 
analyze if we can distinguish between different kinds of SFHs 
with the exceptional photometric coverage that the combi- 
nation of GALEX and ALHAMBRA data provides. 

We build another set of BC03 templates considering 
that the SFH varies exponentially with time, with values 
t sfh = 0.0001, 0.001, 0.01 0.1, 1, 2 5, 10, and 50 Gyr. In this 
process we adopt again a ISalpeterl l|l955t ) IMF distributing 
stars from 0.1 to 100 Mq and select a constant value of 
Z=0.4Zq for t he metallicity. Dust attenuation is also in- 
cluded with the lCalzetti et all l|2000l l law and parametrized 
through the color excess in the stellar continuum, E S (B~V), 
for which we select values ranging from to 0.7 in steps 
of 0.05. We also include int ergalact i c med ium absorption 
adopting the prescription of iMadaul l)l995T K Figures l9l fTOl 
and QT] show the differences found for the SED-derived age, 
stellar mass, and dust attenuation, respectively, when as- 
suming different SFHs. In those plots, the color code in- 
dicates the ratio of the Xr values between the SED-fitting 
results with different SFHs. In those cases where that ra- 
tio is close to one, the SED-fitting procedure is not able 
to discriminate between different kinds of SFHs. Regarding 
age, Figure [9] indicates that, as expected, there is a degen- 



eracy between age and tsfh, with ages systematically older 
when tsfh decreases. Furthermore, this tendency is more 
accused for galaxies whose SED-derived ages with constant 
SFR are higher. For tsfh \Gyr, the ratio of the Xr values 
obtained with different SFHs are very similar to the unity 
and, therefore, although the SED-derived ages are different 
for different values of tsfh, the SED-fitting procedure is not 
able to distinguish between those SFHs. This produce an un- 
certainty in the age of the galaxies. The ages derived with 
templates associated to constant SFH are typically higher 
than those obtained with templates of finite values of tsfh- 
In this way, the ages derived in the previous section for our 
LBGs should be considered as upper values. 

For tsfh < 1 Gyr, the ratio between the Xr values is 
similar to one for the majority of the galaxies but there is 
a population of them (those for which the age derived with 
constant SFR is higher) whose Xr indicate that the tem- 
plates associated to tsfh < 1 Gyr are much worse fitted. 
Only in these cases we are able to distinguish between dif- 
ferent kinds of SFHs and, thus, as a general trend, we ob- 
tained that the SEDs of our oldest GALEX-selected LBGs 
are better fitted with BC03 templates associated to SFH 
with tsfh \Gyr, including the case of constant SFR. 

Regarding the stellar mass, Figure [TO] indicates that 
templates with lower values of tsfh tend to give lower values 
than those associated to higher values of tsfh- In this way, 
the stellar masses derived with constant SFR and reported in 
the previous section for our GALEX-selected LBGs should 
be understood as an upper limits. Again, as for the age, 
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Figure 11. Differences between the SED-derived dust attenuation when performing SED fittings with BC03 templates associated to 
constant SFR and BC03 templates associated to exponentially declining SFH with finite values of the SFH time scale. The color of the 
points in each figure are related to the ratio in the \r values between the SED-fitting results when considering the different kinds of 
SFHs. The values of such Xr ratios corresponding to each color are indicated by the color bars. Values close to one indicate that the 
templates associated to different SFHs fit the photometry with the same accuracy and, therefore, it is not possible to distinguish between 
different kinds of SFHs. Each plot is associated to one value of the SFH time scale, as indicated in each vertical axis. 



even with the good photometric coverage of GALEX and 
ALHAMBRA we are not able to distinguish between the 
different SFH scenarios in most cases. In Figure [10] we rep- 
resent a horizontal line in each plot which corresponds to a 
deviation of 0.3 dex with respect to the case of constant SFR. 
This is the maximum median difference that we find and it is 
a measurement of the uncertainty of the SED-derived stellar 
mass. Figure [TT] shows that there is also a degeneracy be- 
tween the SED-derived dust attenuation and the SFH time- 
scale. In this case, the tendency between dust attenuation 
and tsfh is not as clear as those for age and stellar mass. 
The differences between the results with constant SFR and 
finite values of tsfh are within a factor of two, although for 
most galaxies both estimations agree. 

Summarizing, even with the good photometric coverage 
of the combination of GALEX and ALHAMBRA data we 
are not able to distinguish between different SFHs for most 
of the cases. We can only say that older galaxies (in terms of 
the SED-derived age with templates associated to constant 
SFR) are better fitted with templates associated tsfh 1? 
lGyr. The differences between the SED-derived properties 
when assuming different kinds of SHFs are a measurement 
of the SED-derived procedure. 



5.3 SFR-stellar mass plane 

One important parameter in the study of physical prop- 
erties of galaxies and their evolution is the specific star 
formation rate (sSFR). This quantity is defined as the ra- 



tio between the star formation rate and the stellar mass, 
sSFR=SFR/M„, and it is a measurement of the present over 
past star formation activity. A tight relation between the 
star formation rate and the stellar ma ss has been found 
at different redsh i fts in several works dSalim et al.l 
Elbaz et all 120071; iNoeske et ail 120071; iDaddi et a! 



20101 ; iKarim et aljfeoill ; (Elbaz et alJ ' l2oij ). This has en 



2007; 



2007 



Pannella et al.ll2009l; IRodighiero et ajj|201fj iGonzalez et all 



abled the definition of a 'main sequence' (MS) of galaxies 
defined by a specific value of the sSFR. Galaxies in the 
MS are opposite to the idea of 'starburst galaxies', which 
are those galaxies whose nature makes them deviate from 
the main sequence towards higher values of the sSFR. The 
characteristic value of the sSFR for the main sequence of 
galaxies has been reported to change with redshift (see for 
example lElbaz et al.l (|201lT > for a recent discussion of this 
evolution with deep Herschel FIR data). Previous studies 
show that the sSFR increases with increasing redshifts at all 
masses and that the sSFR of massive galaxie s is lower than 
that for less massive galaxies at any redshift dFeulner et alj 
| 2005l;lErb et alfeOOfj; Dunne et al.ll2009l ; lDamen et al.ll200a 
IRodighiero et alj|2010r~ Despite the number of studies an- 
alyzing the relation between stellar mass and SFR, there 
is still some controversy, mostly regarding the slope of the 
sSFR-M* relation. 

In the left plot of Figure [12] we show the dust-corrected 
total SFR for our GALEX-selected LBG against their stel- 
lar mass. Along with the data poi nts, the 'ma i n sequ ence' 
of galaxies at z ~ 1 taken from lElbaz et al.l l|201ll ) and 
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Figure 12. Left ■plot: Dust-corrected total SFR for our GALEX-selected LBGs as a function of their stellar mass. Red symbols are 
related to the dust corrected total SFRs obtained with the SED-derived value of dust attenuation (as explained in Scction l3.1H . Orange 
symb ols are relat e d to t he dust corrected total SFRs obtained with the dust attenuation derived from the UV continuum slope using 
the iMeurer et alj Jl999T) law. Dark grey and clear grey straight lines are the linear fits to the red a nd orange symbols, respectively. 
For comparison, we also show the SFR- mass relations for the main sequence of galaxies reported in lElbaz et alj l|201ll ) (light green 
continuous line) and Noe ske et alj i2007fi (dark green continuous line) at a simila r redshift range th an that for our GALEX-selected 
LBGs. The dashed light red line represent twice the value of the main sequence of lElbaz et alj l|201lT l. Right plot: Specific SFR of our 
GALEX-selected LBGs as a function of their stellar mass. The meaning of the colors and the fits is the same that those in the left plot. 



iNoeske et all (|2007l ) are also represented. In that plot, red 
symbols indicate the dust-corrected total SFR obtained fol- 
lowing the procedure explained in Section 13.11 Some other 
works, recover the dust attenuation of galaxies employing 
the UV continuum slope and a relation bet ween the UV 
contin uum slope and dust attenuation, as the IMeurer et all 
(jl999h law. The values of dust attenuation obtained in this 
way can be also used for obtaining the dust-corrected total 
SFR. The values found following this procedure are repre- 
sented in Figure [12] with orange symbols. It can be seen 
that whatever the dust-correction method applied, there is 
a clear relation between SFR and mass in the sense that 
more massive galaxies have also higher SFRs. 



Right plot of Figure [12] shows the relation between the 
sSFR and the stellar mass of our GALEX-selected LBGs. 
Again, we represent the points associated to the dust- 
corrected total SFR obtained from the SED-derived dust 
attenuation a nd also those derived from the UV continuum 
slope and the IMeurer et all |l999) law. It can be seen that 
galaxies with larger masses tend to have lower values of 
the sSFR, independently of the dust-correction method em- 
ployed. By fitting a linear relation to the points we obtained 
that the sSFR-mass relation of our GALEX-selected LBGs 
has a slope of -0.49 ± 0.03 and -0.38 ± 0.02 when dust- 
corrected with the SED-derived dust attenuation and the 
UV continuum slope, respectively. 



6 SIZES AND MORPHOLOGY 

In this section we study the morphology and the physical 
sizes of our GALEX-selected LBGs. To this aim, we use I- 
band ACS images taken from the HST observations of the 
COSMOS field, that for our redshift range correspond to a 
rest-frame band ranging from approximately 3500 to 4250 
A. Out of the whole sample of LBGs, 897 have ACS in- 
formation. We download 8"x8" ACS cut-outs centered at 
the position of each source from the cut-outs service avail- 
able on the NASA/IPAC Infrared Science Archive 0. In or- 
der to study the morphology of the galaxies we carry out 
a visual inspection of each cut-out and classify them into 
six groups: Disk-galaxies, compact galaxies, chain galaxies 
(CH), clump cluster (CC), interacting/merging galaxies and 
Irre gular galaxies. For the definition of CH and CC we fol- 
low [EimeaiinelaD (|2009h and consider a galaxy as irreg- 
ular when its morphology does not fit in any of the other 
groups. As a result of the visual classification, we find that 
the majority (69%) of our GALEX-selected LBGs are disk- 
type galaxies, 11% have indications of interactions or merg- 
ing, 7% are irregular, 6% are compact galaxies and a mi- 
nority are CH or CC galaxies, with a fraction of 5% and 
2%, respectively. It should be noted that a visual classi- 
fication is always a very subjective procedure, and there- 
fore, the previous percentages are approximated. However, 
what it is quite clear is that the dominant morphology of 
our GALEX-sele cted LBGs at z ~ 1 is the disk-like class. 
IWolf et all (|2005l ) studied the contribution to the UV lumi- 

2 http:/ /irsa. ipac.caltech.edu/data/COSMOS/index_cutouts. html 
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Figure 13. Distribution of the effective radii (left plot) and the 
Sersic indices (right plot) for our GALEX-selected LBGs with 
available ACS images. Histograms have been normalized to their 
maxima in order to clarify the representations. 



nosity density of different morphologies of galaxies at z ~ 0.7 
by combining hi gh-resolution images from the GEMS survey 
jRix et alj|2004) with redshift s and SEDs from COMBO-17 
survey (|Wolf et al.ll200ll . 120041 ). They report that seemingly 
normal disk galaxies are the largest contributors to the UV 
luminosity density. This is compatible with the result that 
we find in the present work, since LBGs are among the 
brightest UV galaxies and most of them have disk-like mor- 
phologies. Furthermore, the domination of disk-like galax- 
ies morpholo gy in LBGs at z ~ 1 is in agreement with 
the result of iBurgarella et al.l (|200fJ ) for LBGs at a simi- 
lar redshift range than ours, who found that 75% of LBGs 
are compatible with such m orphology. However, ou r result 
contrasts with that found in iBasu-Zvch et all l|2011i) . where 
most LBGs appear to have irregular morphologies and only 
a few have disk-like morphologies. This difference is likely 
due again (see Section [S]) to the inclusion in their sample of 
galaxies at redshift higher than ours. 

In order to obtain the physical sizes of our GALEX- 
selected LBGs we use the previous 8"x8" I-band ACS 
cutouts and carry out fit s to their radial light curves with 
GALFIT ( Peng et al. |2010T ). In this step, we consider Sersic 
profiles (|sersidll96&r ). which can be described as: 



S(r) = S e exp 



r 

ReS 



l/r, 



- 1 



(3) 



where E e is the pixel surface brightness at the effective ra- 
dius i? e ff and n is the concentration parameter or Sersic 
index. The effective radius is the radius which encloses half 
the light of the galaxy. To make this definition true, th e de- 
pendent variable k is coupled to n l|Peng et alj|2010h . For 
each input galaxy, GALFIT provides the effective radius (in 
pixels) and the Sersic index. In order to convert the effective 
radius in pixels into the physical size in kpc we employ the 
ACS pixel scale and the assumed cosmology for calculating 
the " / pix at the redshift of each galaxy. Figure [13] shows 
the distribution of the effective radius and the Sersic index 
for our GALEX-selected LBGs. The median effective radius 
for our LBGs is 2.48 kpc. The values of the Sersic indices 
for LBGs are compatible with most of them being disk-like 
galaxies. This is in agreement with the results of the visual 
morphological analysis. 

Shown in Figure [14] are the relations between the phys- 
ical sizes of our GALEX-selected LBGs and their dust- 



corrected total SFR and stellar mass. It can be seen that 
there is a tendency between the physical size and both SFR 
and stellar mass: larger galaxies tend to form stars faster, 
and have higher stellar masses. The size-stellar mass rela- 
tion has also reported to occur in L B Gs a t a similar and 
higher redshift ranges. iMosleh et all (|2012T l find that the 
stellar mass-size for LBGs persists up to z ~ 5. 



7 COLOR-MAGNITUDE DIAGRAM 

An important tool to analyze the properties of our GALEX- 
selected LBGs is their location in a color-magnitude dia- 
gram. Traditionally, this kind of diagram has been used to 
separate local galaxies between non-SF galaxies earlier than 
the Sa morphological type and SF galaxies later than Sb 
in morphological type. In a color space, the former tend to 
populate a red sequ ence and the lat t er are located in th e 
so-called blue cloud ijHogg et al.ll2002l : IStrateva et alfcOOlh . 
This behavior translates into a bimodal distribution of the 
color of galaxies which allows either to study the nature 
of different samples of galaxies by looking at their posi- 
tion in the color space and to look for galaxies with dif- 
ferent SF nature, by imposing conditions to their location 
in such a diagram. Furthermore, this bimodality in the lo- 
cal universe has been p roven to apply a t higher redshifts , 

1. 12 



at least up to z ~ 1 .6 ilBell et alil2004l; IMcoI et al.l [201ll: 
Williams et al.l 120091: iFranzetti et al.l 120071: ICirasuolo et al. 



2007; iTavlor et all 12003 : IWeiner et alj|2005l : iBlanton et al 
2003). Following this idea, we plot in Figure [T31 the location 
of our GALEX-selected LBGs in a color-magnitude diagram 
associated to the magnitudes in the u and r broad-band fil- 
ters of the SDSS survey. The u and r and absolute r mag- 
nitudes are obtained by taking the templates which fit their 
observed SEDs best and convolving them with the transmis- 
sion of the u and r SDSS filters shifted in wavelength accord- 
ing the redshift of each source. Along with those points we 
also repre sent a sample of SDSS local galaxies taken from 
the DR7 l|Abazauan et all I2009T ) and a sample of z-phot- 
selected galaxies at z ~ 1 taken from the ALHAMBRA sur- 
vey. The sample of local galaxies comprises all the galaxies in 
the SDSS whose spectroscopic redshifts are below 0.035. In 
this case, the magnitudes plotted are those that we extract 
from the photometric catalogs of the SDSS survey. At such 
a low redshifts there is no need of K-correction. In order to 
build the sample of galaxies at 2 ~ 1 we select all the galax- 
ies in the ALHAMBRA survey (in all the already observed 
fields) whose photometric redshifts are around that value 
and whose observed Ks-band magnitudes ar e brighter than 
22 ma g, similar to the limi ts employed in IWilliams et al.l 
( 2009) or lTavlor et~al"1 (|2009T ). By using the optical and near- 
IR photometry of the ALHAMBRA survey we fit their SEDs 
with BC03 templates and obtain their u — r colors, and ab- 
solute r-band magnitude in the same way than for LBGs. It 
can be seen in the left plot of Figure [15] that the bimodal- 
ity which is seen in the local universe is also clearly present 
at z ~ 1. This result also indicates the power of the AL- 
HAMBRA survey in characterizing the CMD of galaxies at 
different redshifts. 

The majority of our GALEX-selected LBGs are located 
in the blue-cloud of galaxies at their redshift, indicating that 
these kinds of galaxies are blue and active SF galaxies, as 
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Figure 14. Left: Effective radius against the dust-corrected total SFR for our GALEX-sclcctcd LBGs with available ACS images. Grey 
squares represent the median value of the dust-corrected total SFR for each considered bin of effective radius. Grey straight line represents 
a linear fit to the grey squares. Right: Effective radius against the SED-derived stellar mass for our GALEX-selected LBGs (red dots) 
with available ACS images. Grey squares represent the median value of the SED-derived stellar mass for each considered bin of effective 
radius. Grey straight line represents a linear fits to the grey squares. 



expected from their selection criteria in the UV. A minor- 
ity of LBGs and are shifted toward the red sequence or are 
located between the blue cloud and the red sequence, the 
so-called green valley. This position does not necessarily in- 
dicates that these galaxies are non-SF. Actually, it can also 
occur that these galaxies have redder optical colors because 
of either a significant amount of dust which is attenuating 
their bluest emitted light, and/or there is an important con- 
tribution of old stellar populations in their SEDs, being more 
evolved systems. To clarify this issue we show in the right 
plot of Figure [Pol the position of our GALEX-selected LBGs 
in the CMD as a function of SED-derived dust attenuation 
and age. We arbitrarily consider two subclasses within the 
LBGs: Those with age > 1200 Myr (old-LBGs) and those 
with age < 1200 Myr and E„(B - V) > 0.4 (dusty-LBGs). 
It can be clearly derived from Figure [15] that those LBGs 
which are located over the green valley or near the red se- 
quence are old-LBGs and dusty-LBGs, whereas those LBGs 
whose ages are bellow 1200 Myr and have low/intermediate 
(E S (B — V) < 0.4) dust attenuation are located over the 
blue cloud. 



8 COMPARISON TO HIGH-REDSHIFT LBGS 

In this section we analyze the differences/similarities be- 
tween LBGs at z ~ 1 and z ~ 3 in order to study whether 
the Lyman break selection criterion selects different kind of 
galaxies at different redshifts. As commented in Section [421 
high-redshift LBGs tend to be intrinsically brighter than 
those studied in this work as a consequence of the usage 
of magnitude-limited observations. If we want to compare 
galaxies which are selected at different redshift with a simi- 
lar selection criterion we have to limit the luminosities of the 



samples to the same range. LBGs are selected through rest- 
frame UV colors and, therefore, in this section we only con- 
sider those LBGs at z ~ 1 whose rest-frame UV luminosities 
are comparable to those for high-redshift LBGs. This sam- 
ple was defined in Section [4.21 as UV- bright LBGs and it is 
formed by 65 galaxies. According to their SEDs, UV-bright 
LBGs are less dusty, have higher SFR, are more massive, 
and have bluer UV continuum slope than the whole popu- 
lation of LBGs at z ~ 1. However, there is no significant 
difference in their ages, which are mostly below 500 Myr for 
both populations. 

Figure [16] shows the distribution of the SED-derived 
age, dust attenuation, and stellar mass for our UV-bright 
LBGs (orange histograms) and for a sample of LBGs at 
z ~ 3 studied in lPapovich et a l. (2001) (green histograms). 
UV-bright LBGs at z ~ 1 have ages mainly distributed 
between 1 and 400 Myr with a median value of 171 Myr, 
whereas LBGs at z ~ 3 are younger galaxies with a median 
age of 36 Myr. Furthermore, as it can be seen in the his- 
togram shown in Figure I16l at z ~ 1 there is a presence of 
older stellar population (with ages mainly between 150 Myr 
and 400 Myr) than at z ~ 3. This could indicate that the 
galaxies selected with the Lyman break criterion at z ~ 1 are 
in a later evolutionary stage than those at z ~ 3. It should 
be noted that, as pointed out in many previous works and in 
Section [5] the uncertainties of the SED-derived stellar age 
are usually high and, therefore, the age evolution found can 
be the consequence of diverse factors: a) using different SFHs 
in the analysis of the SED of the galaxies, b) employing pho- 
tometric information with different wavelength coverage, c) 
the degeneracy between dust attenuation, metallicity, and 
age, etc. Therefore, the previous age evolution of the age of 
LBGs with redshift should be confirmed with further stud- 
ies. 
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Figure 15. Left: Locus of our GALEX-selected LBGs (red dots) in a color-magnitude diagram. For comparison, we represent with blue 
dots the location in such a diagram of a sample of local galaxies selected from SDSS survey. Furthermore, we represent with black contours 
the typical CMD for galaxies at z ~ 1 obtained from a general population of galaxies at that redshift taken from the ALHAMBRA 
survey. These contours show the location of the blue cloud, green valley, and red sequence at that redshift and help in the discussions 
shown in the text. Right: As in the left plot, but we segregate between old-LBGs (brown) and dusty-LBGs (orange). Old-LBGs are those 
LBGs whose ages are above 1200 Myr while dusty-LBGs are those LBGs whose dust attenuation are higher than E S (B — V) = 0.4. 



The dust attenuation distribution of LBGs at z ~ 1 
seems to contain lower values than the distribution at high 
redshift, although both have median values of E S (B — V) = 
0.25. The typical uncertainty of the SED-derived dust at- 
tenuation in our work is AE 3 (B — V) = 0.1 (see Section 
[5]). This value along with the similar median values of the 
distributions at z ~ 1 and z ~ 3 prevent us to constraining 
any evolution of the median values of dust attenuation of 
LBGs with redshifts. This can be an effect of the procedure 
employed. A SED fitting technique is not precise enough 
to constrain an evolution of dust attenuation with redshift 
and others techniques should be employed. Direct measure- 
ments of dust emission of LBGs in the FIR could give clues 
for addressing this issue. 

Regarding stellar masses, it can be seen that the distri- 
butions at z ~ 1 and z ~ 3 span within a similar range. 
The median values of the stellar mass of our UV-bright 
LBGs and LBGs at z ~ 3 are log (M»/M ) = fO.O and 
log(M*/M©) = 9.7, respectively. This difference is similar 
to the typical uncertainties of the stellar mass determina- 
tions in this work and, therefore, we cannot constrain an 
evolution in the stellar mass of LBGs with redshift either. 
The median value of the stellar m ass found in the pre sent 
work is between to those reported in|Magdis et al.1 (|201Gh for 
lRAC-8/im detected and IRAC-8/im faint LBGs at z ~ 3, 
(log M*/M Q ) =11, and (log M»/M©} =9, respectively. 

Regarding the UV continuum slope, UV-bright 
LBGs have a median value of /3 = —1.44. This 
value is l arger (redder) than tho s e reported at higher 



IBouwens et al.1 l2011al : ICastellano et al] l2012f ). indicating 
that LBGs at lower redshifts tend to be redder in the UV 
continuum than those at higher redshifts. The UV contin- 
uum slope is not a parameter associated directly to the BC03 
templates, but it is obtained once the best-fitted BC03 tem- 
plate for each galaxy is known. Consequently, this parameter 
is more insensitive to the different SFH adopted for building 
the BC03 template and, therefore, is a good and accurate 
indicator of the evolution of LBG with redshift. 



redshifts iLehnert fc Bremerl [20031: IBouwens et ah 
IHathi et alJl200St IBouwens et al.ll2009l : IWilkins et al.1 



2006 



2011 



IMosleh et al.1 (|201 il l study the redshift evolution of the 
physical sizes of samples of LBGs and other UV/sub-mm- 
selected galaxies at different redshifts and find tha t their 
size increases with decreasing redshift. IMosleh et af] (|2012T ) 
study the redshift evolution of LBGs from z ~ 1 up to z ~ 7 
and find that the median size of LBGs at a given stellar 
mass increases toward lower redshifts. T he median size of our 
sample of UV-bright LBGs is 2.92 kpc. IMosleh et al.l (|2012T ) 
study the redshift evolution of LBGs considering galaxies 
in two bins of stellar mass, 8.6 < log (M*/Mq) < 9.5 
and 9.5 < log (M*/M Q ) < 10.4. As indicated before, the 
median value of the stellar mass of our UV-bright LBGs 
is log (M,/Mq) — 10.0. The median value of the size of 
UV-bright LBGs is sligh tly lower than that presented in 
IMosleh et all (|2012T i and IMosleh et all (|201lT) for the cor- 
responding stellar mass range. This small difference is likely 
due that, although they also work with GALEX-selected 
LBGs, they consider galaxies located within 0.6 < z < 1.4, 
whereas we limit the redshift of our sample to z > 0.8. The 
inclusion of galaxies at lower redshifts might increase the 
median value of the size, explaining the difference found be- 
tween the two works. 
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Figure 16. Distribution of the SED-derived age, dust atten- 
uation, and stellar mass for our UV-bright LBGs (orange his- 
tograms) and h i gh-re dshift LBGs (green histograms) taken from 
IPapovich et ail l|200ll l. Histo grams have been normalized to their 
maxima in order to clarify the representations. 



It is important to remark again that there is a funda- 
mental difference in the selection criteria of LBGs at differ- 
ent redshifts (see Section ??). Given that the combination 
of GALEX and ALHAMBRA (including IRAC for a sub- 
sample) provides very accurate photometric redshifts of our 
galaxies, we do not need any extra condition to rule out in- 
terlopers. However, this is not the case at high redshifts. At 
z > 2 the photometric redshifts of the sources are not accu- 
rate enough to ensure a proper cleaning of the sample from 
interlopers. As a consequence, additional criteria should be 
applied. These extra criteria usually involve limits in the 
observed op tical colors of the sa mples. For example, as it is 
discussed in lMadau et~aH (|l996t ) for their sample of F300W 
dropouts, the application of extra criteria rules out interlop- 
ers which are at another redshift range, but they also missed 
galaxies at the proper redshift at the same time. The missing 
galaxies tend to be redder than those included in the final 
sample, either because they are older or are attenuated by 
dust. Therefore, it is clear than the additional criteria em- 
ployed at high redshift get rid of certain kinds of subclasses 
of galaxies at ea c h red shift. In contrast, in our wor k, as in 
IBurgarella et al.l |2006l ) and iBurgarella et al.l (|2007m , we do 
not apply any extra observed optical color criteria and there- 
fore we include in the sample all kind of galaxies which have 
a break between FUV and NUV filters, regardless their age, 
dust attenuation or optical colors. 

All previous results seem to indicate that LBGs at z ~ 1 
tend to be older, bluer in their UV continuum, and have 
larger sizes than those at higher redshifts. Therefore, SED 
fitting and morphological studies indicate that LBGs at 
lower redshifts are in a later evolutionary stage than those 
at higher redshifts. It should be considered that these differ- 
ences between LBGs at different redshifts have been found 
by comparing the results obtained in this work with other 



from previous studies done by other authors. Therefore, this 
comparison might suffer from uncertainties coming from the 
different methods employed in each work or slightly different 
selection criteria. For example, the use of BC03 templates as- 
sociated to different SFHs or different photometric coverage 
of the SED of galaxies at different redshifts. Then, in order 
to properly characterize the evolution of LBGs with redshift 
a more precise work should be done where the photometric 
SEDs of the galaxies are sampled over the same rest-frame 
wavelength range and the selection criteria and the proce- 
dures employed for the analysis of the physical properties of 
the galaxies are as similar as possible. 



9 CONCLUSIONS 

In this work we have analyzed the SED-fitting derived phys- 
ical properties of a sample of 1225 GALEX-selected LBGs 
at 0.8 < z < 1.2 by using a combination of UV and 
optical/near-IR data coming from GALEX observations and 
the ALHAMBRA survey, respectively. ALHAMBRA uses a 
set of 20 medium-band (width ~ 300Jl) optical and the three 
classical near-IR JHKs filters for covering the observed opti- 
cal SED of galaxies in an unprecedented way. This provides 
a good sampling of both the UV continuum slope and the 
Balmer break, increasing significantly the accuracy of the 
results of the SED-fitting technique. We defined LBGs as 
those galaxies that have a difference of color greater than 
1.5 magnitudes between the FUV and NUV filters of the 
GALEX satellite. Our main conclusions are as follows: 

(i) According to SED fitting with BC03 templates built 
assuming a constant SFR, Salpeter IMF, and metallicity 
Z = OAZq, GALEX-selected LBGs at z ~ 1 are young 
galaxies with ages mostly below 300 Myr, with median dust 
attenuation of E S (B — V) of 0.20, and a median UV contin- 
uum slope of -1.53. As a consequence of the selection criteria 
used they are UV-bright objects with UV-uncorrected SFR 
of about 2.0 Mq /yr. When dust-correcting their rest-frame 
UV luminosity, their total SFR turns out to have a median 
value of 46.4 M©/yr. Combining the total SFRs and ages, 
we find that GALEX-selected LBGs have a median stellar 
mass of log (M*/M ) = 9.74. Only 2% of the galaxies selec- 
tion with the Lyman break selection criterion have an AGN 
according to their X-ray emission. 

(ii) LBGs with higher stellar masses have higher to- 
tal SFRs and lower values of the specific SFR. The anti- 
correlation between the specific SFR and stellar mass sup- 
ports the downsizing scenario, where more massive galaxies 
have formed their stars earlier and faster than galaxies with 
lower stellar mass. 

(iii) Morphologically, LBGs at z ~ 1 are mostly disk-like 
galaxies (about 69%), while the remaining are interacting, 
compact or irregular systems in much lower percentages. 
This is confirmed by their Sersic indices, which are typi- 
cally below 0.5. The median effective radius for our GALEX- 
selected LBGs at z ~ 1 is 2.48 kpc. Bigger galaxies tend to 
have higher total SFR and stellar masses. 

(iv) In a color-magnitude diagram, most GALEX-selected 
LBGs are located over the blue cloud at their redshift, which 
indicates that they are active SF galaxies. Some LBGs are 
located over the green valley or near the red sequence. They 
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turn out to be the dustiest and/or oldest galaxies in the 
sample, as signs of more evolved systems. 

(v) Comparing with their high-redshift analogs, we find 
that the galaxies selected through the Lyman break crite- 
rion at z ~ 1 seem to be in a later evolutionary stage than 
those at high-redshift. LBGs at lower redshifts are bigger, 
have more contribution of older stellar population to their 
SEDs, and are redder in their UV continuum than their 
high-redshift analogues. However, we do not find any signif- 
icant difference in the distributions of stellar mass or dust 
attenuation for LBGs at high and intermediate redshift. 
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